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Executive  Summary 


The  proposed  research  brought  forth  the  novel  idea  of  constructing  intersubband 
lasers  on  Si-based  material  systems,  which  could  potentially  offer  the  great  opportunity 
for  monolithic  integration  of  Si-based  photonic  devices  with  advanced  Si  electronics.  The 
material  systems  under  investigation  include  quantuam  structures  of  SiGe/Si  for  mid- 
to  far-infrared  laser  applications  and  of  Si/ZnS  for  the  near-infrared. 

Specifically,  we  have  studied  the  scattering  of  heavy-holes  by  optical  phonons  in 
SiGe/Si  quantum  well  structures,  which  dominates  the  lifetimes  of  heavy-hole  subbands 
in  SiGe/Si  intersubband  lasers.  Intrasubband  and  intersubband  scattering  rates  of  heavy 
holes  are  obtained  due  to  confined  nonpolar  optical  phonons  in  a  Sii_xGex  quantum 
well  with  Si  barriers.  Guided  and  interface  Ge-Si  and  Ge-Ge  modes  and  unconfined 
Si-Si  modes  are  considered.  A  continuum  model  is  used  for  the  two  components  of  the 
ionic  displacement  of  confined  vibrations:  the  uncoupled  s-polarized  TO  mode  and  the 
hybrid  of  the  LO  and  p-polarized  TO  modes.  The  guided  mode  is  obtained  using  the 
model  of  a  quantum  well  with  infinitely  rigid  barriers,  and  the  interface  mode  is  derived 
from  the  hydrodynamic  boundary  conditions.  While  the  total  intersubband  scattering 
rates  are  reduced  as  a  result  of  confinement,  the  opposite  is  found  for  the  intrasubband 
scattering.  Depending  on  the  well  width  and  Ge  content,  the  intersubband  scattering 
rates  are  reduced  by  a  factor  of  two  to  four  with  respect  to  their  values  for  no  confinement. 
Thus,  one  would  expect  comparable  enhancement  in  the  intersubband  lifetimes  crucial 
to  the  population  inversion  in  a  Sii_zGez/Si  intersubband  laser. 

We  have  also  calculated  the  valence  subband  engeries  and  wavefunctions  in  Si/ZnS 
quantum  well  structures.  We  have  also  established  models  for  the  confined  optical 
phonons  in  Si/ZnS  superlattices.  Optical  phonon  modes  in  Si  and  ZnS  layers  are  totally 
confined  within  their  respective  layers  since  both  layers  can  be  treated  as  infinitely  rigid 
with  respect  to  the  other  layer.  Since  there  are  no  associated  electric  fields  with  nonpolar 
optical  phonons  in  Si  layers,  only  a  mechanical  boundary  condition  needs  to  be  satisfied 
for  these  nonpolar  optical  modes  at  the  Si-ZnS  interface.  The  optical  phonons  in  Si 
layers  can  be  described  by  guided  modes  consisting  of  an  uncoupled  s-TO  mode  and  a 
hybrid  of  LO  and  p-TO  modes  with  no  interface  modes.  In  ZnS  layers,  a  continuum 
model  hybridizing  the  LO,  TO  and  IP  modes  is  necessary  to  satisfy  both  the  mechanical 
and  electrostatic  boundary  conditions  at  the  heterointerface.  A  numerical  procedure  is 
provided  to  determine  the  common  frequency  between  LO,  TO,  and  IP  modes.  This  is 
a  new  procedure  for  obtaining  the  eigen  modes  of  a  mixed  polar-nonpolar  heterosystem. 
Analytical  expressions  are  obtained  for  the  ionic  displacement  and  associated  electric 
field  as  well  as  scalar  and  vector  potentials.  The  established  model  can  be  used  to  calcu¬ 
late  the  phonon  scattering  process  of  carriers  in  Si/ZnS  heterostructures  and  ultimately 
provide  a  theoretical  basis  for  studying  the  Si/ZnS  intersubband  lasers. 
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Part  I.  A  study  of  confined  nonpolar  optical  phonons  and  their  interaction 
with  carriers  in  SiGe/Si  quantum  well  structures 

I.  INTRODUCTION 

The  possibility  of  lasing  due  to  intersubband  transitions  in  quantum  well  (QW) 
structures  is  crucially  dependent  on  the  lifetimes  of  the  involved  subbands.  In  an  ear¬ 
lier  communication [1]  lifetimes  were  calculated  for  a  Sii_xGes/Si  multiple  quantum  well 
(MQW)  structure  due  to  acoustic  and  nonpolar  optical  phonon  scattering.  It  was  pointed 
out  that,  because  of  the  absence  of  polar  optical  scattering  for  silicon-based  systems,  the 
lifetimes  and  their  differences  are  consistently  an  order  of  magnitude  larger  than  for  the 
GaAs/AlGaAs  system,  and  do  not  show  the  marked  decrease  when  the  intersubband 
energy  exceeds  the  optical  phonon  energy.  For  this  case,  the  scattering  was  calculated 
for  propagating  phonons  corresponding  to  the  average  composition,  assuming  no  con¬ 
finement.  As  pointed  out  in  that  paper,  it  is  expected  that  this  treatment  of  phonon 
scattering  leads  to  an  overestimate  of  the  intersubband  transition  rates,  and  therefore 
provides  a  conservative  approach  in  determining  the  subband  lifetimes.  In  order  to  more 
accurately  determine  the  subband  lifetimes,  it  is  necessary  to  take  into  account  confine¬ 
ment  effects  on  optical  phonons  in  heterostructures.  In  a  Sii_xGex/Si  QW  structure, 
one  can  expect  Si-Si,  Ge-Si,  and  Ge-Ge  optical  modes  in  the  alloy.[2]  The  Ge-Si  and 
Ge-Ge  modes  tend  to  be  confined  by  the  Si  barriers  as  has  been  seen  in  Raman  scat¬ 
tering  experiments.  [3]  The  phonon  confinement  gives  rise  to  guided  and  interface  modes 
which  will  scatter  the  carrier  less  for  intersubband  transitions,  thus  resulting  in  revised 
subband  lifetimes.  This  is  due  to  the  discrete  spectrum  of  the  guided  modes  and  the 
weak  intersubband  scattering  of  interface  modes.  However,  the  intrasubband  scattering 
process  tends  to  be  enhanced  due  to  confinement. 

There  has  been  a  great  deal  of  effort  in  dealing  with  the  issue  of  phonon  confine¬ 
ment  in  heterostructures. [2]-[21]  Most  of  published  literature  have  focused  on  confined 
polar  optical  phonons,  because  a  large  body  of  the  heterostructures  are  constructed  from 
polar  semiconductor  materials.  Recently,  the  demonstration  of  an  infrared  intersubband 
transition  laser  in  the  InGaAs/AlGaAs  material  system[22]  has  renewed  interest  in  in¬ 
vestigating  the  possibility  of  constructing  lasers  within  the  Ge/Si  material  system,  which 
eventually  would  allow  monolithic  integration  of  optical  components  with  advanced  Si 
microelectronics. 

A  self-consistent  continuum  theory  of  confined  nonpolar  optical  phonons  has  been 
developed  for  an  infinite  plate  with  free  boundary  conditions.  [4]  In  the  present  paper, 
we  extend  that  continuum  theory  to  examine  the  confined  nonpolar  optical  phonons 
of  Ge-Si  and  Ge-Ge  vibration  modes  with  proper  mechanical  boundary  conditions  in  a 
Sii-xGex/Si  QW.  Furthermore,  we  use  these  results  to  estimate  the  heavy-hole  scattering 
rates  by  the  nonpolar  optical  phonons.  To  the  best  of  our  knowledge,  there  has  not  been 
any  work  in  estimating  the  scattering  rates  due  to  the  nonpolar  optical  phonons  in 
Sii_sGex/Si  heterostructures  taking  into  account  phonon  confinement. 
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The  three  different  vibration  modes,  namely,  Si-Si,  Ge-Si,  and  Ge-Ge,  are  consid¬ 
ered  separately  and  are  given  proper  weights  in  heavy-hole  scattering  calculation.  The 
Si-Si  mode  propagates  freely  throughout  the  structure,  while  the  Ge-Si  and  Ge-Ge  modes 
are  confined  by  the  barriers,  resulting  in  guided  and  interface  modes.  Approximations 
made  in  the  course  of  the  calculations  for  the  guided  and  interface  modes  will  be  dis¬ 
cussed.  We  will  examine  both  intrasubband  and  intersuband  scattering  processes.  We 
will  also  compare  these  results  with  those  assuming  no  confinement.  We  will  show  that 
indeed  the  latter  assumption  overestimates  the  scattering  rates  for  intersubband  tran¬ 
sitions,  but  tends  to  underestimate  the  intrasubband  process.  The  difference  is  largely 
attributed  to  the  heavy-hole  interaction  with  interface  modes  which  can  be  neglected 
for  the  intersubband  process,  but  contributes  significantly  to  the  intrasubband  process. 
The  scattering  rates  due  to  various  scattering  modes  will  be  investigated  as  a  function 
of  the  well  width  and  as  a  function  of  the  alloy  composition  of  the  well. 

II.  PHONON  CONFINEMENT 

Phonon  confinement  occurs  due  to  the  lack  of  overlap  of  the  bulk  frequency 
dispersions  in  the  adjacent  heterostructure  materials.  Thus  in  a  Si/Ge/Si  short  period 
superlattice,  the  higher  frequency  optical  modes  of  Si  correspond  to  an  frequency  gap 
in  the  Ge  layer  and  are  therefore  confined.  The  Ge  modes  in  the  Ge  layers  are  strictly 
resonant  with  the  Si  acoustic  modes  and  are  ‘quasi-confined’,  but  the  displacement 
pattern  is  similar  to  that  of  a  true  confined  mode  and  indeed  the  Raman  intensities 
are  comparable. [2]  The  Raman  spectrum  has  revealed  that  there  exist  Ge-Ge,  Ge-Si, 
and  Si-Si  vibration  modes  in  the  Sii_iGez  alloy  layer. [3]  The  Ge-Ge  vibrations  have 
phonon  energy  of  ZlAmeV  and  the  Ge-Si  vibrations  of  b0.8meV.  These  two  modes 
of  the  Sii_zGex  well  are  confined  by  the  Si-Si  modes  of  the  barrier  (phonon  energy 
64.3meF);  however,  the  Si-Si  optical  modes  of  the  well  are  unconfined,  as  are  the  Si-Si 
optical  modes  of  the  barrier. 

One  of  the  simplest  conceptual  model  is  to  treat  the  QW  system  with  infinitely 
rigid  barrier.  The  boundary  condition  to  be  satisfied  is  the  vanishing  of  the  ionic  dis¬ 
placement  of  all  confined  vibration  modes.  This  is  an  assumption  of  strict  confinement, 
yielding  guided  modes  of  the  confined  phonons.  This  assumption  needs  to  be  partially 
relaxed  in  order  to  admit  interface  modes.  As  pointed  out  in  the  continuum  theory,  [4] 
the  ionic  displacement  of  confined  vibrations  has  two  components:  one  is  the  hybrid  of 
the  LO  and  p-polarized  TO  (p-TO)  modes,  and  the  other  is  the  uncoupled  s-polarized 
TO  (s-TO)  mode.  These  modes  are  defined  as  follows:  if  we  consider  a  (x,  z)  plane 
containing  the  normal  to  the  layers  and  the  phonon  wavevector  Q,  then  . 

Q  =  qxex  +  qzez  (1) 

where  ex  and  ez  are  unit  vectors.  The  p-TO  mode  has  its  displacements  normal  to  Q 
and  in  the  (x,  z )  plane,  while  the  s-TO  displacements  are  normal  to  Q  and  perpendicular 
to  the  (x,  z)  plane  (||ey).  A  description  of  the  s-TO  mode  is 

uy  =  ei9lI(Aei9‘z  +  Be-iq‘z ),  (2). 
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(3) 


while  the  hybrid  of  the  LO  and  p-TO  modes  is  given  by 

Ux  =  eiQxX[qx(CeiqLZ  +  De~iqLZ)  +  qT{EeiqTZ  +  Fe~iqTZ)\, 

Uz  =  eiQxX[ qL(CeigLZ  -  De~iqLZ)  -  qx(EeiqTZ  -  Fe~iqTZ) ]. 

The  2-components  of  the  LO  and  TO  wavevector  have  been  distinguished  by  qi  and  qr, 
respectively. 

The  above  choice  for  the  ionic  displacement  guarantees  that 

V  x  ut  =  0  and  V  •  Ut  =  0,  (4) 

which  hold  for  isotropic  materials  (assumed  here),  allowing  u  to  be  decomposed  into  LO 
and  TO  components, 

u  =  uL  +  ur.  (5) 

Since  the  LO  and  TO  modes  must  have  the  same  frequency  to  be  effectively 
coupled,  they  have  to  satisfy,  according  to  the  bulk  LO  and  TO  dispersions, 

uj2  =  uJq  —  P2L{ql  +  q\)  =  Wq  -  Pt(qI  +  <?r)>  (6) 

where  (3l  and  fir  are  the  velocities  of  LO  and  TO  dispersions,  respectively. 

Guided  Modes 

In  the  case  of  strict  confinement,  the  boundary  condition  of  course  is  that  the 
displacements  u  vanish  at  the  boundaries  z  =  (0,  L ).  For  the  s-TO  mode,  this  leads  to 

TVK 

uy  =  AetqxX  sin(qzz),  with  qz  —  —  (7) 

Lj 

where  n  —  1,2,--  *.  This  mode  does  not  mix  with  other  modes,  nor  does  it  give  rise  to 
the  interface  mode  discussed  later.  The  s-TO  mode  pattern  for  qz  =  7 r/L  is  shown  in 
Fig.  1(a)  with  the  well  width  L  —  40 A. 

For  the  coupled  LO  and  p-TO  modes,  applying  the  boundary  conditions  of  ux  =  0 
and  uz  =  0  at  z  —  (0,  L)  to  Eq.(3)  leads  to  a  system  of  four  constant-coefficient  linear 
equations  for  C,D,E,  and  F.  The  vanishing  of  the  determinant  of  the  4x4  constant 
coefficient  matrix  gives  rise  to  the  following  relation, 

{qrqL  +  ql)2  sin (qTL)  sin (qLL)  =  2qTqLql[cos{qT  -  ?l)T  -  1]-  (8) 


Eq.(8)  leads  to  solutions  of  guided  modes  consisting  of  coupled  phase-matched 
LO  and  TO  modes  with  wavevectors 


TlL'K  . 

qL  =  —f~  and  qr 


TVilT 

IT 


(9) 
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where  til  =  1, 2,  •  •  •,  nr  =  3, 4,  •  •  •,  and  nr  —  ni  =  2, 4, 6,  •  •  This  choice  of  quantum 
numbers  is  due  to  the  constraint  Eq.(6)  since  fir  <Pl- 


Figure  1.  Vibration  patterns  for  (a)  the  guided  s-TO  mode,  (b)  the  ‘sine’  solution,  (c) 
the  ‘cosine’  solution  of  the  guided  p-TO  and  LO  modes,  and  (d)  the  coupled  p-TO  and 
LO  interface  mode  for  the  well  width  of  40 A. 


Two  sets  of  guided  mode  solutions  emerge  and  they  are  given  in  either  sine  or 
cosine  form  referring  to  the  z-component  of  the  ionic  displacement.  The  ‘sine’  solution 
is 


ux  =  2CetqxXqx[cos(qLz)  -  cos(gTz)], 

a 2 

uz  =  2  iCeiqxX[qL  sin(qLz)  +  —  sin(qTz)], 

Qt 


(10) 


and  a  ‘cosine’  solution 


ux  =  2 iCeiqxX[qx  sin(qLz)  +  sin(gTz)], 

Qx 

uz  =  2CetqxXqL[cos(qLz)  —  cos(qtz)]. 


(11) 


The  hybrid  patterns  of  the  p-TO  and  LO  modes  with  qi  =  7r/ L,  q x  =  Zn/L  are  shown  in 
Fig.l(b,c)  for  the  ‘sine’  and  ‘cosine’  solutions,  respectively,  with  the  well  width  L  =  40A. 
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It  is  worthwhile  to  point  out  the  similarities  between  the  allowed  values  of  qz,  qi, 
qT  and  qx  for  the  guided  modes  in  the  case  of  infinitely  rigid  barriers  assumed  here  and 
that  of  an  infinite  plate  with  free  boundary  conditions.  [4]  But  the  displacement  patterns 
differ  dramatically  between  the  two  situations. 

As  pointed  out  by  Mori  and  Ando[l3],  Nash[14],  and  Register[15],  a  sum  rule  can 
be  derived  for  polar  optical  scattering  relating  the  form  factors  of  the  confined  modes  to 
that  of  the  ‘bulk’  modes.  Such  a  sum  rule,  as  pointed  out  by  Nash[14],  arises  from  the 
completeness  and  orthogonality  of  the  vibration  sets.  As  derived,  this  applies  to  polar 
optical  scattering  accompanied  by  a  macroscopic  electrical  potential.  As  is  well  known, 
the  latter  removes  the  degeneracy  of  the  zone  center  optical  phonons  and  the  carrier 
interacts  only  with  the  LO  mode.  However,  for  the  nonpolar  optical  phonons  of  interest 
in  the  present  work,  both  TO  and  LO  modes  interact  with  the  carrier.  As  discussed 
earlier,  the  s-TO  mode  does  not  mix  with  the  other  modes  and  can  be  described  in  terms 
of  an  orthogonal  and  complete  set  of  functions  given  by  Eq.(7).  However,  the  p-TO  and 
LO  modes  are  coupled  as  given  by  Eqs.(10)  and  (11),  with  the  constraint  given  by  Eq.(6) 
that  the  frequencies  be  equal,  taking  account  of  the  different  frequency  dispersions.  The 
LO  and  p-TO  cannot  be  treated  independently  because,  having  chosen  their  forms  to 
satisfy  Eq.(4),  the  vanishing  boundary  boundary  conditions  at  z  =  0  and  z  =  L  cannot 
be  satisfied  for  u l  and  ut  independently.  As  a  consequence  of  the  coupling,  in  satisfying 
the  boundary  conditions  for  the  guided  modes,  not  all  values  of  wavevectors  ( <7x><7z,>9r ) 
are  allowed  and  so  the  set  is  not  complete.  Specifically,  the  allowed  values  of  (qL,  qr) 
require  that  qx  be  discrete,  not  continuous,  and  be  given  by 

2  _  /7T>2  pTnT  ~  PWl 

where  ni  and  nr  have  been  defined  in  Eq.(9).  Further,  two  modes  which  differ  in  only 
one  of  (qi„  qr)  are  not  orthogonal.  In  summary,  in  comparison  with  the  guided  modes  for 
polar  optical  phonons,  even  though  the  guided  modes  constitute  a  complete  orthogonal 
set  for  the  s-TO  modes  in  the  y-direction,  the  linear  combination  of  p-TO  and  LO  modes 
in  the  ( x ,  z)  plane  are  neither  complete  nor  orthogonal.  Hence,  a  sum  rule  cannot  be 
applied  to  the  current  basis  set  for  nonpolar  optical  phonons. 

Interface  Modes 

The  assumption  of  strict  confinement,  that  the  vibration  amplitudes  are  zero 
at  the  interfaces,  rules  out  the  possibility  of  interface  modes.  In  general  such  modes 
exist  and  are  of  increasing  importance  in  carrier  scattering  intrasubband  transitions  in 
comparison  with  guided  modes.  The  assumption  of  strict  confinement  must  therefore  be 
relaxed  to  admit  such  modes.  It  is  obvious  that  the  amplitudes  of  interface  modes  at  the 
boundaries  should  remain  small  in  order  to  keep  the  assumption  of  strict  confinement 
yielding  the  guided  modes.  In  fact,  we  need  only  to  relax  the  vanishing  assumption  of 
the  x-component  of  the  displacement  to  admit  interface  modes  while  maintaining  the 
zero  boundary  condition  on  the  z-component.  A  similar  boundary  condition  has  been 
employed  in  treating  the  polar  optical  phonons  for  a  single  GaAs/ALAs  QW.[10] 
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Applying  the  hydrodynamic  boundary  conditions[5]  at  the  interfaces  to  the  LO 
and  TO  modes  independently  and  considering  the  fact  that  the  optical  phonon  frequency 
in  the  barrier  region  is  greater  than  that  in  the  well  region,  we  obtain  only  a  solution 
with  even  z-component  for  the  LO  mode, 


r 


Ul  =  < 


v. 


CeiqxX(qxex  -  qL2ez)e  iqL2Z 

CeiqxX(— )1/2(  — — qxex  -  qL2ez) 
.  P2  qiPi  t 
—CetqxX(qxex  +  qL2ez)e'q^z~V 


z  <  0, 

0  <  z  <  L, 

z  >  L, 


(13) 


where  the  subscripts  1  and  2  refer  to  well  (Sii_xGez)  and  barrier  (Si)  regions,  respectively, 
Pi  is  the  material  density,  and 

rii  =  w?/u>2  -  1,  (14) 

where  cuf  is  the  T-point  optical  phonon  frequency  in  layer  i  (=  1,2).  For  both  Ge-Si  and 
Ge-Ge  confined  modes  we  have  u2  >  uii.  The  LO  wavevector  in  the  well  region  qLl  =  0 
and  that  in  the  barrier  region, 

The  LO  mode  solution  is  of  two-dimensional  bulk  type  with  constant  amplitudes  in  the 
well  region  propagating  with  wavevectors  parallel  to  the  interfaces. 


Both  odd  and  even  solutions  emerge  for  the  TO  mode,  but  the  boundary  condition 
that  the  z-component  of  ur,-f  Ur  vanishes  at  the  interfaces  admits  only  the  even  solution, 


f  DeiqxX[ctr2ex  +  qxez)e  iqT2Z 
Pi  -  • 1 


u T-  < 


where  the  wavevectors  for  the  TO  mode,  according  to  Eq.(6) 


Qti 

=  (Pk-  i)a2 

{P2  )Qx 

(well), 

QT2 

=  #  - 1  )<&  +  #«L 

PT2  PT2 

(barrier). 

z  <  0, 


DeiqxX  (-)1/2 - 

Kp2)  cos(9T1L/2) 
x  [-iqn  sin  qTi(z  -  L/2)ex  +  qx  cos qTi(z  -  L/2)ezJ  0  <  z  <  L, 
[  -DeiqxX{qT2ex  -  qxez)eiq^z~V 


z  >  L, 


(16) 


(17) 


Combining  Eqs.(13)  and  (16)  and  applying  that  the  z-component  of  ux,  +  Ur  vanishes 
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at  the  interfaces,  lead  to  D  =  ( quijq^C ,  and  finally  the  interface  mode  is  obtained, 


u  = 


Ceiq*x[{qxe 


’'IMS*  — 


.  Qx  . 

+  (-<ZL2e-‘*“*  +  gi2e-,9T2Z)ez] 

— —  sin  qn  (z  -  L/2)]ex 
Pi'  " m  pi  “  cos(oTiL/2j  9l 
cosgn(z  —  I>/2)  _ 

•  r  /  rA  1  ZJ 


1  V  ^ 
CeWx*(*)l/2{[ 
Oo 


z  < 


+  cos(qTiL/2) 

CeiqxX{[-qxeiqL^z~L)  -  q~T2  eiqT2(z~L)]ex 

+  [-9£2ei9“(z_L)  +X9z,2ei9r2(z_i:)]e2} 


0  <  2  <  L, 


z  >  L. 


(18) 


Fig.  1(d)  shows  the  pattern  of  an  interface  mode  in  the  well  region  responsible 
for  the  scattering  process  from  subband  2  to  1  by  emitting  a  Ge-Si  optical  phonon  with 
L  =  40 A.  It  can  be  seen  from  Fig.  1(d)  that  the  n-component  of  the  displacement  is 
approximately  constant  within  the  well  region  and  the  z-component  vanishes  at  both 
interfaces. 


The  above  form  of  the  interface  mode  derived  from  the  hydrodynamic  boundary 
conditions  is  different  from  that  of  the  surface  mode  for  the  infinite  free  plate[4]  and  that 
for  the  interface  mode  for  polar  materials.  The  difference  from  the  former  case  arises 
from  the  different  boundary  conditions  since  the  free  surface  requires  that  the  dilation 
stress  perpendicular  to  the  surface  and  the  shear  stress  across  the  surface  vanish.  The 
difference  from  the  latter  is  primarily  due  to  the  fact  that  the  zone-center  optical  phonon 
frequency  in  the  barriers  is  greater  than  that  in  the  well,  which  is  contrary  to  the  situation 
in  GaAs/AlGaAs  material  systems. 


The  dispersion  of  bulk  Si-Si  mode  has  little  overlap  with  that  of  Ge-Si  mode. 
The  scattering  process  involves  only  long  wavelength  in-plane  wavevectors,  and  the  z- 
components  of  phonon  wavevectors  in  the  barriers,  qi2  and  qn,  are  approximately  the 
dimension  of  the  Brillouin  zone.  This  is  because  the  dispersion  relations  overlap  only  at 
large  phonon  wavevectors  in  the  barrier  region.  Since  qn  =  0  and  qx  «  0,  the  interface 
mode  frequency  ui  «  u>i,  therefore,  772/771  >>  1.  The  major  scattering  contribution  from 
the  interface  mode  is  from  the  displacement  in  the  well  region  where  the  heavy-holes 
are  confined.  The  amplitude  of  the  LO-component  is  much  greater  than  that  of  the 
TO-component  in  Eq.(18)  for  0  <  z  <  L,  since 


V2P2  gL2gn 

Vi  Pi  9* 


(19) 


considering  qn  also  small.  Therefore,  the  exact  values  of  qn  and  qr 2  are  not  crucial 
in  determining  the  scattering  rates  because  the  TO  mode  contribution  is  at  least  two 
orders  of  magnitude  less  than  the  LO  mode.  The  above  displacement  pattern  Eq.(18) 
has  been  used  in  evaluating  the  interface  mode  scattering  for  the  confined  Ge-Si  and 
Ge-Ge  modes. 


9 


III.  HEAVY  HOLE  SCATTERING  RATES 


is  [23] 


The  nonpolar  optical  phonon  interaction  Hamiltonian  involving  a  heavy  hole 


Mi  +  M2 


(20) 


where  Mi  and  M2  are  the  masses  of  the  two  atoms  in  the  unit  cell,  and  D  is  the  optical 
deformation  potential. 


The  normalization  of  the  displacement  amplitudes  is  carried  out  through  the 
approach  of  equating  the  energy  of  the  vibration  mode  with  that  of  a  simple  harmonic 
oscillator[9]  as 

*2=§ru'-u<i2’  <2i> 

where  L  is  the  well  width,  S  is  the  sample  surface  area  (in  (x,  y)  plane),  Q  is  the  volume 
of  the  unit  cell,  and  x  is  the  normal  coordinator  of  the  oscillator. 
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Figure  2.  Heavy  hole  band  structure  and  scattering  processes. 

The  heavy-hole  band  offset  is  calculated  taking  into  account  the  compressive 
strain  in  the  Sii_IGeI  well  region. [24]  The  heavy-hole  bands  are  decoupled  from  the 
light-hole  and  split-off  bands  at  k  =  0  and  can  be  treated  independently.  [25]  The  heavy- 
hole  energy  levels  and  envelope  wavefunctions  are  obtained  by  the  finite  square  well 
model  as  shown  in  Fig.2.  The  heavy-hole  state  can  be  characterized  by  |k,  n  >  with  the 
in-plane  momentum  k  and  subband  index  n.  In  the  approximation  of  constant  effective 
mass  for  heavy  holes,  the  matrix  element  for  the  transition  from  state  |k,  n  >  to  Ik',  v!  > 


10 


due  to  nonpolar  optical  phonon  scattering  is 


<  k',n'\H\k,n  >=  { 


(S'T0) 

Ili[n{u0)  +  WtJS5i 

2piu0SLAc(qL, 9t)  k  ±qx,k 

•  [DXG*.  (qL,  qT)  +  DzGznnl  {qL,  qT)]  (hybrid) , 


(22) 


for  the  s-TO  mode  and  the  hybrid  of  the  LO  and  p-TO  mode,  respectively.  n(u0)  is 
the  number  of  optical  phonons  at  thermal  equilibrium,  and  the  upper  and  lower  signs 
refer  to  phonon  absorption  and  emission,  respectively.  The  three  components  of  the 
optical  deformation  potential,  Dx,  Dy,  and  Dz  are  assumed  equal  to  D0  —  D/y/Z  in  the 
calculation,  in  view  of  the  assumption  of  isotropy.  The  Kronecker  symbol  indicates  the 
in-plane  (x,  y )  momentum  conservation.  The  normalization  factors  are  given  by 

^■a{Qz)  —  y  f  U  U  dz  (s-TO), 

LJ\  ,L  (23) 

Ac{qL,  qr )  =  j  Jq  (Kux  +  u*zuz)dz  (hybrid). 

The  Gnn< -functions  contain  envelope  wavefunctions,  tpn  and  ip'n,  from  which  interference 
effect  can  be  obtained.  Specifically, 


for  the  s-TO  mode,  and 


Gln'(qz)  =  f  Mn'Uydz, 
J  0 


fL 

Gxnn,{qL,qT)=  ipni>n>uxdz, 

J0L 

GZnn.(qL,qT)  =  /  Tpn^n'Uzdz, 
Jo 


(24) 


(25) 


for  the  hybrid  of  LO  and  p-TO  mode.  It  should  be  noted  that  given  the  proper  dis¬ 
placement  expressions,  Eq.(22)  for  the  matrix  element  is  valid  for  both  the  guided  and 
interface  modes  discussed  above. 


Obviously,  depending  on  the  alloy  composition  in  the  well,  the  three  modes  (Si- 
Si,  Ge-Si,  and  Ge-Ge)  will  have  different  interaction  strengths  with  the  carriers  in  the 
QW  structure.  Specifically,  each  interaction  should  be  weaker  than  that  for  the  case 
where  it  is  the  only  existing  mode  and  therefore  needs  to  be  accounted  for  properly.  A 
crude  model  to  approximate  the  relative  strength  of  each  individual  mode  is  to  assign 
a  weight  in  the  calculation  of  the  phonon  scattering  rate.  If  we  assume  bonds  formed 
between  Si-Si,  Ge-Si,  and  Ge-Ge  are  purely  random,  then  we  can  give  the  following 
weights  according  to  the  Ge  content,  x,  to  each  of  the  vibration  modes  in  the  Sii^Ge* 
well:  wsi  =  (1  —  x)2  for  Si-Si  mode,  waesi  =  2x(l  —  x)  for  Ge-Si  mode,  and  wae  =  x 2 
for  Ge-Ge  mode. 


ll 


Applying  the  Fermi  golden  rule,  we  obtain  the  scattering  rate  due  to  the  guided 


modes, 


wfK,  [nK)-fl/2Tim.  1  G^l 

2h2Pluj0L  ^  Aa 

^•m^[n(a;0)  +  l/2Tl/2]D0V  |G*W,  +  G^,]2 

2  h2Plu0L  ^qL'gT  Ac 


(s-TO) 

(hybrid), 


(26) 


where  we  have  assumed  that  for  the  intersubband  process  (n  ^  n')  the  heavy  holes  are 
scattered  from  the  bottom  of  their  original  subbands,  and  for  the  intrasubband  process 
( n  =  n ')  the  heavy  holes  have  just  enough  kinetic  energy  (equal  to  the  Si-Si  optical 
phonon  energy)  to  emit  an  optical  phonon  to  reach  to  the  bottom  of  the  same  subband. 
Wj  ( j  =Ge-Si,  Ge-Ge)  is  the  weight  assigned  to  a  particular  vibration  mode  and  m*hh  is 
the  heavy-hole  effective  mass. 


For  the  s-TO  guided  mode,  the  summation  in  Eq.(26)  is  over  all  qz  =  mrjL 
limited  by 

£  <  (f  ?  -  £,  (27) 

ai 

where  Oj  is  the  lattice  constant  and  qx  is  given  by  the  constraint  of  the  in-plane  momen¬ 
tum  conservation.  For  the  hybrid  guided  mode,  the  summation  is  restricted  to  those 
combinations  of  qi,  —  nin/L  and  qr  =  nr^/L  which,  according  to  Eq.(6),  yield  discrete 
values  of  qx  satisfying  the  in-plane  momentum  conservation.  It  should  be  pointed  out 
that  an  unique  final  state  heavy-hole  wavevector  k'  (initial  k  =  0)  will  result  if  the 
phonon  frequency  dispersion  is  neglected,  corresponding  to  a  constant  optical  phonon 
energy.  Then  the  in-plane  momentum  conservation  would  be  impossible  since  the  dis¬ 
crete  values  of  qx  will  not  be  able  to  exactly  match  the  unique  value  of  k'.  However, 
taking  into  account  of  the  phonon  dispersion,  a  region  of  heavy-hole  wavevector  k'  will 
be  obtained  and  the  scattering  process  is  permitted  as  long  as  qx  lies  in  that  region. 


There  is  no  s-TO  interface  mode,  and  the  scattering  rate  due  to  the  hybrid 
interface  mode  can  be  given,  similarly, 


+  1/2  T  l/2]r>g 

piUgLAc 


l  G;„.+G'nf, 


(28) 


where  G*n,  and  Gznn,  are  functions  of  qri,  Ql2>  and  qr 2  with  qn  =  0.  The  determination 
of  the  phonon  wavevectors  in  the  well  and  baxrier  regions  have  been  discussed  in  section 
II.  Eqs.(26)  and  (28)  are  used  to  evaluate  both  the  intrasubband  and  intersubband 
transitions  due  to  the  confined  phonons  of  Ge-Si  and  Ge-Ge  modes.  However,  since  the 
Si-Si  vibration  is  unconfined,  the  calculation  of  its  scattering  rate  is  trivial.  [26] 

Intrasubband  Scattering 


The  intrasubband  scattering  rate  was  calculated  assuming  that  the  transition 
originated  from  the  heavy-hole  state  with  a  kinetic  energy  equal  to  the  Si  phonon  energy 
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to  the  bottom  of  the  same  subband.  Therefore,  only  the  process  of  phonon  emission 
is  considered  in  the  calculation.  All  results  are  obtained  at  room  temperature.  The 
structure  parameters  are  varied  within  the  limits  for  producing  metastable  strained 
Sii-sGe*  alloy  on  the  Si  substrate. [27]  Fig.3  shows  the  intrasubband  transition  rates 
(1-1)  within  the  ground-state  heavy-hole  subband  (n  =  1)  due  to  the  Ge-Si  and  Ge-Ge 
guided  and  interface  modes  as  a  function  of  the  well  width  for  a  Sio.sGeo.s/Si  QW.  It 
can  be  seen  that  for  the  intrasubband  process  the  strength  of  interface  mode  scattering 
is  about  same  order  of  magnitude  as  that  of  corresponding  guided  mode.  This  is  in 
contrast  with  the  intersubband  process  discussed  later.  The  scattering  rates  of  both 
guided  and  interface  modes  for  Ge-Si  and  Ge-Ge  vibrations  increases  with  the  increase 
of  the  well  width  when  it  is  narrow  (L  <  25 A).  This  is  due  to  the  increase  of  the 
interference  G-function  with  n  =  n'  =  1  given  in  Eqs.(24)  and  (25)  as  the  envelope 
function  for  subband  1  becomes  more  confined  in  the  well  region  for  larger  well  width. 
However,  the  G-function  increase  becomes  negligible  as  the  well  width  increases  further 
since  the  envelope  function  has  mostly  been  confined  within  well  region.  Examining 
the  interface  scattering  rate  given  in  Eq.(28),  one  would  expect  it  to  decrease  with  the 
further  increase  of  the  well  width,  which  is  clearly  demonstrated  in  Fig.3  for  both  Ge-Si 
and  Ge-Ge  phonon  interactions.  A  similar  decrease  of  the  interface  scattering  rate  in  a 
GaAs/Ala;Gai_*As  QW  has  been  reported. [7] 


Figure  3.  Intrasubband  scattering  rates  in  the  lowest  subband  (n  =  1)  due  to  Ge-Si  and 
Ge-Ge  guided  and  interface  modes  as  a  function  of  well  width  with  x  =  0.5. 

For  the  guided  mode,  the  scattering  rate  is  given  by  summing  the  contributions 
from  the  s-TO  and  hybrid  of  LO  and  p-TO  modes.  However,  the  hybrid  interaction 
with  the  heavy  holes  is  much  weaker  than  the  s-TO  mode,  because  of  the  requirements 
on  the  hybrid  mode  to  simultaneously  satisfy  the  in-plane  momentum  conservation  and 
the  frequency  dispersion  Eq.(6).  This  excludes  the  interaction  of  most  hybrid  modes 
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with  the  heavy  holes,  leaving  the  s-TO  modes  as  the  dominant  scattering  mechanism. 
In  fact,  the  number  of  allowed  s-TO  vibration  modes  will  increase  with  the  increase  of 
the  well  width.  -  It  is  easy  to  see  from  Eq.(24)  that  for  the  intrasubband  process  (n  =  n') 
the  interference  G- function  is  nonzero  only  for  even  s-TO  modes  (nqt  =  1, 3, 5,  •  •  •), 
and  the  largest  contribution  is  from  the  lowest  order  mode  nqz  =  1  with  higher  order 
modes  being  at  least  two  orders  of  magnitude  less.  The  guided  mode  scattering  for  the 
intrasubband  process  is  therefore  practically  a  single  s-TO  mode  interaction  and  behaves 
similarly  to  the  interface  mode.  The  difference  between  the  Ge-Si  and  Ge-Ge  vibrations 
is  mainly  due  to  the  different  weights  assigned  to  them  according  to  the  Ge  content  in 
the  SisGei-x  alloy.  Our  calculation  also  showed  that  the  (2-2)  intrasubband  scattering 
process  has  very  similar  behavior  as  the  interface  mode  scattering  and  has  a  slightly 
higher  scattering  rate  than  the  guided  mode,  but  overall  smaller  rates  compared  to  the 
(1-1)  process. 

Fig.4(a)  shows  the  total  scattering  rate  as  a  function  of  the  well  width,  which 
is  given  by  summing  the  contributions  from  the  Ge-Si,  Ge-Ge,  and  Si-Si  (unconfined) 
vibrations.  Also  shown  in  Fig.4(a)  is  the  result  calculated  neglecting  the  phonon  con¬ 
finement  for  comparison.  It  can  be  seen  that  the  phonon  confinement  actually  enhances 
the  intrasubband  scattering  process.  However,  the  opposite  is  true  for  the  intersubband 
process  to  be  shown  below.  Fig.4(b)  shows  the  same  result  as  a  function  of  the  Ge 
content  in  the  alloy  for  the  well  width  of  40 A.  The  difference  increases  with  the  increase 
of  the  Ge  content  since  the  contributions  from  the  confined  phonons  will  increase  while 
that  from  the  unconfined  Si-Si  vibration  decreases. 


Figure  4.  Total  intrasubband  scattering  rates  in  the  lowest  subband  (n  =  1)  due  to  Ge-Si, 
Ge-Ge  confined  modes  and  Si-Si  unconfined  modes,  compared  to  the  bulk  intrasubband 
scattering  rate  assuming  no  confinement,  (a)  as  a  function  of  well  width  with  Ge  content 
x  =  0.5;  (b)  as  a  function  of  x  for  a  well  width  of  40 A. 
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Intersubband  Scattering 


The  intersubband  scattering  rate  was  calculated  assuming  that  the  transition 
originated  from  the  bottom  of  a  subband  n  with  zero  kinetic  energy  to  another  sub¬ 
band  ( n '  /  n).  Both  phonon  absorption  and  emission  processes  are  considered  in  the 
calculation.  Fig.5(a)  shows  the  (2-1)  scattering  rate  with  the  Si-Si,  Ge-Si,  and  Ge-Ge 
components  as  a  function  of  well  width  for  a  Sio.sGeo.s/Si  QW.  Subband  2  appears  at 
the  well  width  of  20A.  The  total  scattering  is  obtained  by  summing  the  contributions 
derived  from  the  Ge-Si,  Ge-Ge,  and  Si-Si  vibrations  including  both  guided  and  interface 
modes.  But  the  interface  mode  scattering  is  at  least  two  orders  of  magnitude  weaker 
than  the  guided  mode  for  the  intersubband  process.  Similar  weakness  of  the  interface 
mode  has  been  shown  in  a  GaAs/AlxGai_j;As  QW.[7]  The  result  assuming  no  confine¬ 
ment  is  also  shown  in  Fig. 5 (a)  for  comparison.  The  sharp  increase  of  the  scattering 
rates  in  the  small  well  width  region  ( L  <  25A)  is  once  again  due  to  the  increasing  con¬ 
finement  in  the  well  region  of  subband  2  which  leads  to  an  increase  of  the  interference 
G-function.  Both  scattering  intensities  from  Si-Si  and  unconfined  phonons  reduce  with 
further  increase  of  the  well  width.  Since  initially  the  narrow  well  width  only  allows  a 
small  number  of  guided  modes,  the  guided  mode  scattering  is  weak  compared  to  the 
Si-Si  scattering  and  the  total  scattering  rate  follows  the  dependence  of  Si-Si  vibration. 

As  the  number  of  allowed  guided  modes  increases  with  increasing  well  width,  the  guided 
mode  scattering  actually  surpasses  the  Si-Si  component,  leading  to  an  increase  of  the 
total  scattering  rate.  The  small  discontinuous  incremental  steps  in  the  Ge-Si,  Ge-Ge, 
and  therefore  total  scattering  curves,  are  due  to  the  discrete  nature  of  the  increase  in 
the  number  of  allowed  guided  modes  as  the  well  width  increases.  The  sudden  drops  in 
the  curves  occur  when  the  energy  separation  between  the  two  subbands  reduces  to  less 
than  one  of  the  three  phonon  energies  corresponding  to  different  vibrations.  Specifi¬ 
cally,  the  drops  occurred  at  L  —  72,  83,  99A  correspond  to  Si-Si  phonon  (64.3meV), 
Ge-Si  phonon  (50.8meV),  and  Ge-Ge  phonon  (37.4meV),  respectively.  It  can  be  seen  in 
Fig.5(a)  that  the  phonon  confinement  reduces  the  scattering  rate  by  a  factor  of  two  to 
four  as  compared  to  the  case  of  no  confinement.  It  is  therefore  important  to  take  into 
account  the  phonon  confinement  effect  in  estimating  the  subband  lifetimes. 

Fig.5(b)  shows  the  (2-1)  scattering  rates  as  a  function  of  Ge  content  in  the 
Sii-xGea,  alloy  for  a  well-  thickness  of  40A.  The  contributions  from  the  guided  modes 
increases  with  the  Ge  content  while  the  unconfined  Si-Si  component  reduces.  As  a  re¬ 
sult,  the  difference  between  the  total  scattering  rate  and  the  rate  obtained  assuming  no 
confinement  widens  with  Ge  content. 

IV.  SUMMARY  AND  DISCUSSION 

In  an  earlier  paper  we  treated  phonon-induced  intersubband  transitions  in  Sii_xGex/Si 
MQW  structures.  [1]  The  absence  of  polar  optical  scattering  resulted  in  subband  lifetimes 
an  order  of  magnitude  larger  than  those  in  MQWs  of  III-V  semiconductors  and  without 
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the  precipitous  decrease  at  the  optical  phonon  threshold.  This  earlier  work  assumed  un¬ 
confined,  three  dimensional  phonons.  Noting  strong  experimental  evidence  for  phonon 
confinement  in  Si1_xGeI/Si  MQWs,  the  present  paper  extends  the  earlier  treatment  to 
the  case  in  which  phonons  are  confined. 


Figure  5.  Total  intersubband  scattering  rates  from  subband  2  to  1  due  to  Ge-Si,  Ge- 
Ge  confined  modes  and  Si-Si  unconfined  modes,  compared  to  the  bulk  intersubband 
scattering  rate  assuming  no  confinement,  (a)  as  a  function  of  well  width  with  Ge  content 
x  =  0.5;  (b)  as  a  function  of  x  for  a  well  width  of  40 A. 

A  single  Sii_xGex  QW  with  Si  barrier  is  considered.  The  decoupled  heavy-hole 
subbands  and  wavefunctions  are  determined  by  the  finite  square  well  model.  Assigning 
statistical  weights  to  the  Si-Si,  Ge-Si,  and  Ge-Ge  vibrations  in  the  QW  assuming  a  purely 
random  bond  model,  the  Ge-Si  and  Ge-Ge  modes  are  treated  as  confined  while  the  Si-Si 
modes  are  extended.  The  confined  modes  are  of  two  kinds,  guided  modes  and  interface 
modes.  The  guided  modes  consist  of  the  s-TO  and  coupled  LO  and  p-TO  modes,  with 
the  boundary  condition  that  the  displacements  vanish  at  the  interfaces.  In  comparison 
with  the  guided  modes  for  polar  optical  phonons,  the  current  choice  of  basis  set  which 
describes  the  coupling  between  the  p-TO  and  LO  modes  are  neither  orthogonal  nor 
complete.  Hence,  the  sum  rule  found  for  polar  LO  phonons  cannot  be  applied  to  the  case 
of  nonpolar  optical  phonons.  The  interface  mode  is  obtained  applying  hydrodynamic 
boundary  conditions  to  the  LO  and  TO  modes  separately  and  then  requiring  that  the 
z-component  of  their  sum  (not  the  x-component)  vanish  at  the  interfaces.  The  LO 
mode  is  a  two-dimensional  bulk-type  solution  with  constant  amplitude  in  the  well  and 
propagating  parallel  to  the  interface;  the  TO  mode  is  an  even  function  (referring  to  the 
x-component)  with  a  sinusoidal  spatial  dependence  so  that  the  boundary  conditions  at 


16 


the  interfaces  can  be  satisfied. 

The  Ge-Ge  optical  phonon  spectrum  in  the  well  region  has  no  overlap  with  that 
of  Si-Si  optical  phonons  in  the  barriers,  but  instead  overlaps  with  that  of  the  acoustic  Si 
phonons.  The  model  for  the  interface  mode  derived  from  the  optical  phonon  dispersion 
relations  using  the  hydrodynamic  boundary  conditions  is  not  without  criticism.  Further 
investigation  of  this  subject  is  underway  to  provide  a  more  accurate  estimate  of  the 
intrasubband  scattering  rates.  But  we  would  still  expect  the  intersubband  process  to 
be  dominated  by  the  guided-mode  scattering,  which  ultimately  determines  the  subband 
lifetimes. 

The  intrasubband  and  intersubband  scattering  rates  are  calculated  as  a  function 
of  the  QW  structure  parameters:  Ge  content,  x,  in  the  Sii_xGex  well  and  the  well  width, 
L.  These  parameters  are  varied  within  the  limits  reportedly  to  produce  the  metastable 
strained  layers  of  Sii_xGex  alloy  on  Si  substrates  without  introducing  a  large  number  of 
misfit  defects.  We  choose  the  Ge  content  of  x  =  0.5  as  we  vary  the  well  width  to  give  a 
fair  weight  of  the  Ge-Si  and  Ge-Ge  confined  modes.  We  use  the  fixed  well  width  of  40A  to 
allow  at  least  two  confined  heavy-hole  subbands  as  we  vary  the  Ge  content.  Our  results 
are  as  follows.  For  intrasubband  scattering,  the  scattering  by  the  interface  mode  is  nearly 
the  same  as  that  for  the  guided  mode  (Fig.  3),  and  the  total  intrasubband  scattering  with 
confinement  is  larger  than  that  without  confinement  (Fig. 4).  However  for  intersubband 
scattering,  the  interface  mode  scattering  is  at  least  two  orders  of  magnitude  smaller  than 
that  for  the  guided  mode.  Finally,  as  can  be  seen  from  Fig.5,  the  total  intersubband 
scattering  rate  with  confinement  is  a  factor  of  two  to  four  smaller  than  that  which  would 
be  obtained  assuming  no  confinement.  Thus,  one  would  expect  a  factor  of  two  to  four 
enhancement  in  the  intersubband  lifetimes  which  determine  population  inversion  in  an 
intersubband  laser. 

Finally,  we  note  that  the  present  investigation  has  considered  the  scattering  of 
heavy  holes  whose  dispersion  is  parabolic.  The  extension  to  light  holes  and  the  split-off 
band,  which  are  coupled  due  to  spin-orbit  coupling  and  strain [25]  and  whose  dispersions 
are  therefore  nonparabolic,  would  require  expension  of  the  theory,  but  with  similar 
results  anticipated. 
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Part  II.  A  study  of  confined  optical  phonons  in  Si/ZnS  superlattices 

I.  INTRODUCTION 

There  is  an  interest  in  moving  the  lasing  from  the  far-  and  mid-infrared  range  to 
the  near-infrared  optical  communications  wavelengths,  A  =  1.3  or  1.55/im[29].  Since  the 
latter  wavelength  corresponds  to  a  photon  energy  of  800meV,  the  Si^Ges/Si  heterosys¬ 
tem  is  inadequate  because  a  maximum  practical  valence  band  offset  of  only  ~  500meV 
can  be  obtained  for  pseudomorphic  Sii-xGej  layers  at  x  =  0.5  ~  0.6.  Therefore,  al¬ 
ternative  large-  bandgap,  nearly  lattice-matched  barrier  materials  for  Si  quantum  wells 
must  be  sought;  materials  with  sufficiently  large  band  offsets  with  respect  to  silicon. 
Possible  candidates  include,  ZnS,  BeSeTe,  CaF2,  SiOs,  Si02,  the  Si/Si02  superlattice, 
and  7-AI2O3,  among  others[29]-[31].  The  Si/ZnS  heterosystem  has  received  the  most 
attention  as  current  advances  in  epitaxy  technology  have  allowed  the  growth  of  het¬ 
erostructures  consisting  of  polar  and  nonpolar  materials[32,  33].  The  lattice  mismatch 
of  cubic  ZnS  with  respect  to  Si  is  only  0.3%.  The  valence  band  offset  has  been  predicted 
theoretically [34]- [37],  while  recent  experiments[38]  show  that  the  value  is  close  to  1.5eV, 
sufficiently  large  to  give  intersubband  energy  differences  in  the  desired  800meV  range. 
MBE  growth  of  ZnS  upon  Si,  and  Si  upon  ZnS  have  been  demonstrated [33],  with  the 
use  of  an  As  monolayer  to  satisfy  the  local  bonding  requirements,  although  the  effect  of 
the  monolayer  on  the  offsets  has  not  been  determined. 

The  possibility  of  population  inversion  and  the  operation  of  intersubband  lasers 
depend  critically  on  the  lifetimes  of  the  involved  subbands.  The  subband  lifetimes  in 
turn  are  determined  by  nonradiative  phonon  scattering  processes.  The  purpose  of  the 
present  paper  is  to  study  the  optical  phonon  modes  and  their  interaction  with  carriers  in 
the  Si/ZnS  system  since  the  optical  phonon  scattering  is  considered  to  be  dominant  in 
the  phonon  scattering  processes.  This  combination  of  materials  is  new,  since  it  consists 
of  both  a  nonpolar  and  polar  semiconductor.  Previous  studies  in  carrier  scattering  by 
confined  optical  phonons  in  heterostructures  have  been  focused  only  on  one  type  of 
phonons,  either  polar[7]-[10]  or  nonpolar[3]-[4].  In  the  current  situation  involving  both 
polar  and  nonpolar  materials,  carrier  scattering  by  both  types  of  phonons  needs  to  be 
considered.  To  the  best  of  our  knowledge,  there  has  not  been  any  reported  investigation 
on  this  mixed  nature  of  optical  phonons,  their  confinement  effect,  and  their  interaction 
with  carriers  in  a  heterostructure.  In  this  paper,  we  will  present  a  theoretical  study 
based  on  the  macroscopic  continuum  model  to  describe  the  confined  optical  phonon 
modes  and  will  use  this  model  to  calculate  the  optical-phonon  scattering  of  heavy  holes 
in  a  heterostructure  consisting  of  polar  (ZnS)  and  nonpolar  materials  (Si),  as  we  are 
interested  in  the  feasibility  of  constructing  an  intersubband  laser  within  the  valence  band 
of  the  Si/ZnS  heterostructure.  This  valence  intersubband  laser  will  likely  be  a  quantum- 
parallel  superlattice  laser[40]  or  a  quantum  cascade  laser.  Our  latest  thinking[41]  is 
that  each  of  the  N  laser  periods  will  consists  of  one  square  Si  quantum  well  containing 
two  active  heavy-hole  subbands.  Due  to  the  non-paxabolicity  of  these  subbands,  we 
believe  that  a  population  inversion,  localized  in  A-space,  can  be  engineered  between 


18 


the  subbands.  In  this  investigation,  we  will  consider  a  simple  superlattice  comprised  of 
alternating  layers  of  Si  and  ZnS,  much  like  the  flat-band  superlattice  of  the  quantum 
parallel  laser  (Fig.2(b)  of  [40]).  The  results  of  this  study  will  provide  the  basis  for  the 
calculation  of  subband  lifetimes  required  to  determine  laser  gain  and  threshold. 

As  described  below  in  greater  detail,  since  the  optical  dispersions  (frequency 
versus  wavevector)  of  the  silicon  (huisi  =  64meV)  and  zinc  sulphide  (hujzns  —  43meV) 
have  no  overlap,  the  optical  phonons  are  assumed  to  be  totally  confined  in  both  materials. 
In  the  silicon  layers,  a  continuum  model  with  double  hybridization  of  the  longitudinal 
optical  (LO)  and  transverse  optical  (TO)  modes  is  used  to  describe  the  vibration  patterns 
of  the  guided  modes  as  described  in  Part  I.  The  only  boundary  condition  that  needs  to 
be  satisfied  in  the  Si  layers  is  the  vanishing  of  the  displacements  at  the  Si-ZnS  interface, 
since  the  ZnS  layers  can  be  considered  as  infinitely  rigid  with  respect  to  the  vibrations 
of  the  Si  layer.  Hence,  there  is  no  interface  mode  in  the  Si  layers.  The  situation  on 
the  ZnS  layers  is  more  complex.  Following,  the  work  by  Ridley[8,  9],  here  a  continuum 
model  is  employed  with  hybridization  of  the  optical  LO,  TO,  and  interface  polariton 
(IP)  modes  needed  to  satisfy  both  the  mechanical  and  electrostatic  boundary  conditions 
at  the  interfaces.  Specifically,  the  electrostatic  boundary  conditions  are  the  continuity  of 
Ex,  the  electric  field  parallel  to  the  interface,  and  the  continuity  of  Dz,  the  displacement 
field  normal  to  the  interface.  The  mechanical  boundary  condition  is  again  the  vanishing 
of  the  optical  displacements  since  the  Si  layers  can  be  considered  as  infinitely  rigid  with 
respect  to  the  vibrations  of  the  ZnS  layers. 

Our  current  work  provides  a  complete  set  of  analytical  expressions  for  the  optical 
phonon  dispersion  relations,  optical  displacements,  and  associated  scalar  and  vector 
potentials.  These  expressions  are  subsequently  used  in  calculating  the  interaction  of 
heavy  holes  with  the  confined  optical  phonons.  In  Section  II,  we  establish  a  continuum 
model  for  the  optical  displacement  modes  in  Si  and  ZnS  layers  satisfying  both  mechanical 
and  electrical  boundary  conditions.  In  section  III,  we  outline  a  numerical  procedure  for 
determining  the  frequency  of  a  ZnS  optical  mode  inducing  the  intersubband  scattering. 
In  Section  IV,  we  describe  the  scalar  and  vector  potentials  associated  with  the  ZnS 
ionic  displacement  modes.  In  Section  V,  we  summarize  and  discuss  our  results  and 
conclusions. 


II.  Mode  Patterns 

A  continuum  model  for  the  optical  modes  in  the  Si/ZnS  superlattice  is  employed. 
Both  mechanical  and  electrical  boundary  conditions  are  satisfied  at  the  heterointerfaces. 
Since  the  optical  dispersion  relations  (frequency  versus  phonon  wavevector)  in  the  two 
bulk  materials  have  no  overlap,  the  phonons  are  taken  to  be  confined  in  their  respective 
materials.  For  the  Si  layers,  the  continuum  model  for  optical  phonons  in  nonpolar 
materials  is  used.  Here  double  hybridization  of  the  LO  (longitudinal  optical)  and  TO 
(transverse  optical)  modes  is  used  to  give  the  vibration  patterns  of  the  guided  modes. 
Since  the  ZnS  layers  are  infinitely  rigid  with  respect  to  the  vibrations  of  the  Si  layers,  only 
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the  mechanical  boundary  condition,  the  vanishing  of  the  displacements  at  the  interfaces, 
has  to  be  satisfied. 

4 

For  the  polar  ZnS  layers,  the  situation  is  more  complex  and  an  alternate  con¬ 
tinuum  model[8,  9]  consisting  of  an  intermixing  of  confined  LO,  TO,  and  IP  (interface 
polariton)  modes  is  needed  in  order  to  satisfy  both  the  electrostatic  and  mechanical 
boundary  conditions.  The  boundary  conditions  which  must  be  satisfied  are  (1)  the  con¬ 
tinuity  of  Ex,  the  component  of  electric  field  parallel  to  the  interface,  (2)  the  continuity 
of  Dz,  the  component  of  the  displacement  vector  normal  to  the  interface,  and  (3)  the 
vanishing  of  the  vector  displacement  u  at  the  interface. 

A.  Modes  in  Si  Layers 


As  discussed  in  detail  in  Part  I,  using  the  boundary  condition  that  u  =  0  at  the 
interfaces  gives  for  the  s-TO  mode 

77,7T 

uy  =  Ae'9xX  sin(qzz),  with  qz  —  —  (29) 

asi 


where  n  =  1,2,***  and  A  is  a  mode  coefficient.  This  mode  does  not  mix  with  other 
modes. 


The  hybrid  LO  and  p-TO  modes  admit  two  classes  of  solutions. 

ux  =  2Bei9xXqx[cos(qLz)  -  cos {qrz)], 
uz  =  2iBeikxX[qL  sin(qLz)  +  —  sin(qTz)] , 


The  ‘sine’  solu- 
(30) 


and  the  ‘cosine’  solution  is 

ux  =  2  iBeiq*x[qx  sin(giz)  +  sin(<jr.z)], 

Qx 

uz  =  2BetQxXqL[cos(qLZ )  -  cos(gr-2)] 


(31) 


where 


Ql  = 


TIlTT 


and  qr  = 


nrTT 


dsi  ^  dsi 

•,  nr  —  til  =  2,4,6,  •  ••, 


(32) 


where  ni  =  1, 2,  •  •  •,  nr  =  3, 4,  •  •  *,  nr  -  nt  =  2, 4, 6,  •  •  •,  and  B  is  a  mode  coefficient. 
No  interface  modes  exist  in  the  Si  layer  because  of  the  boundary  condition  u  =  0. 


The  lowest  s-TO  mode  pattern  in  Eq.(29)  for  qz  =  7T /dsi  is  shown  in  Fig.6(a) 
within  a  Si  layer  of  dsi  =  40A,  while  the  hybrid  patterns  of  the  lowest  p-TO  and  LO 
modes  with  qL  =  n/dsi  and  qr  =  3tt /dsi  are  shown  in  Figs.6(b)  and  6(c)  for  the  ‘sine’  and 
‘cosine’  solutions  given  in  Eqs.(30)  and  (31),  respectively  within  the  same  Si  layer.  The 
strict  confinement  which  requires  the  vanishing  of  ionic  displacements  at  the  boundaries 
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of  Si  layers  is  clearly  demonstrated  for  both  vibration  modes. 


Figure  6.  Vibration  patterns  in  a  Si  layer  with  a  width  of  40 A  for  (a)  the  guided  s-TO 
mode,  (b)  the  ‘sine’  solution,  and  (c)  the  ‘cosine’  solution  of  the  guided  p-TO  and  LO 
modes. 


B.  Modes  in  ZnS  Layers 

.  The  boundary  conditions  are  the  continuity  of  Ex,  Dz,  and  the  vanishing  of  u  at 
the  interfaces.  These  conditions  can  be  satisfied  by  a  unique  linear  combination  of  LO, 
TO,  and  IP  modes  with  common  frequency  and  common  in-plane  wavevector  qx, 


u  —  vllo  +  uto  +  U/p.  (33) 

We  will  use  this  hybrid  expression  to  calculate  the  electrical  interaction  with  carriers 
which  is  considerably  stronger  than  the  optical  deformation  potential  interaction.  We 
need  consider  only  the  displacements  ux  and  uz,  since  uy  associated  with  the  s-TO  mode 
has  no  related  electric  field  and  therefore  does  not  interact  with  carriers  electrically. 
Once  again,  the  expressions  are  obtained  by  taking  the  Bloch  superlattice  wavevector 
Qsl  =  0. 

For  the  LO  mode,  the  ionic  displacements  are 

u*  =  e^‘x-^qx{AL^  +  BLe-^z), 

uz  =  ei(qxX-^qL(ALe^LZ  -  BLe~iqLZ).  ^4' 

which  is  confined  within  the  ZnS  layer  with  a  width  of  dzns >  —dzns/ 2  <  z  <  dzns/% 
The  associated  electric  fields  are 

Ex  =  Po^n  Ez  —  po^zi  (35) 

where 

*  =  (36> 
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with  the  effective  ionic  charge 


e-2  =  MO^loeK—  -  — (37) 

^OO  ^3 

where  M  is  the  reduced  mass,  e0  is  the  permittivity  of  free  space,  Coo,  c*  are  the  high- 
frequency  and  static  permittivities,  and  0  is  the  volume  of  primitive  unit  cell.  The  scalar 
potential  (f>  associated  with  the  electric  field  E  =  — V0  is  in  turn  given  as 

<f>  =  -ip0SqxX-^  {ALe^z  +  BLe~iqLZ) .  (38) 


For  the  TO  mode 

ux  =  ei(9lI_wt)gr(ATe<9rZ  +  BTe~iqTZ), 
uz  =  -e**x-ut>qx{Are?**  -  BTe~iqTZ). 

The  electric  fields  associated  with  this  mode  are  negligible. 

For  the  IP  mode 

ux  =  qp  {Apeiq”z  +  BPe~iq »*), 

uz  =  iei{qxX-ut)qp{Apeiq”z  -  BPe~iq*z) 

The  associated  electric  fields  are 

Ex  =  Pp^xt  Ez  =  Ppiiz) 


(39) 

(40) 

(41) 


where 

Pp  =  Po 

and  lolo  and  u>ro  are  bulk  ZnS  LO  and  TO  optical  phonon  frequencies  at  T  point, 
respectively.  The  electric  fields  associated  with  the  interface  modes  propagate  into  the 
Si  layers  although  they  are  treated  as  infinitely  rigid  and  do  not  contain  ZnS  ionic 
displacement. 

Being  a  transverse  electromagnetic  wave,  there  is  a  vector  potential  A  associated 
with  the  electric  field  E  =  —dA/dt.  Within  the  ZnS  layers, 


,  ,2  ,  ,2 
LO  —  wj>o 

9  2  » 

UlO  -  W TO 


Ax  =  i—ei<'qxX~u>t^qv{APeiqpZ  +  BPe~iq”z), 

lo  (43) 

Az  =  -^ei(qxX~u%(APeiqpZ  -  BPe~iqpZ) 

LO 

While  in  the  Si  layers,  a  similar  expression  can  be  obtained  with  another  set  of  mode 
coefficients,  Ap i  and  Bv\. 

Since  large  in-plane  wavevectors  are  likely  to  be  involved  when  dealing  with  carrier 
transitions  due  to  optical  phonons  between  two  subbands  separated  with  a  relatively 
large  energy,  we  have  endeavored  to  obtain  analytical  dispersions  of  the  LO  and  TO. 
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optical  branches  by  curve-fitting  the  experimental  bulk  phonon  dispersions  for  the  entire 
Brillioun  zone.  The  requirement  for  common  frequency  yields, 


w  =uLO-  0L(ql  +  ql) 
-  wto  —  Pr(<ll  +  Qt) 

_c(<?*  +  <g)1/2 

n(u>) 


(44) 


where  Pl  =  0.808 THz  •  A2  and  fir  —  2.194 THz  •  A2  are  obtained  from  curve-fitting 
the  bulk  ZnS  optical  phonon  dispersions,  c  is  the  velocity  of  light  in  vacuum,  and  n(u>) 
is  the  index  of  refraction.  In  the  above  expressions,  the  frequency  in  the  ZnS  layers 
lies  between  the  ZnS  LO  and  TO  zone  center  frequencies.  Since  loto  <  ^lo,  in  order 
for  the  TO  frequency  to  be  equal  to  a  LO  frequency  qr  must  be  imaginary  qr  =  iq0, 
corresponding  to  a  TO  interface  mode.  The  modes  which  interact  most  strongly  with 
carriers  are  those  with  frequencies  near  the  LO  branch.  For  these  modes,  the  value  of 
q0  is  large,  and  we  can  take  the  approximation 


tanh (q0dzns)  ~  1.  (45) 

In  the  unretarded  limit  (c  oo),  qx  +  g2  «  0  for  the  IP  mode.  Hence,  qp  «  iqx. 

Applying,  at  the  two  interfaces  between  layers  Si  and  ZnS  in  a  period  of  the 
superlattice,  the  conditions  that  ux  and  uz  equal  to  zero  along  with  the  continuity  of  Ex 
and  Dz,  leads  to  eight  simultaneous  equations  involving  the  eight  unknown  mode  coef¬ 
ficients  (Al,  Bl',  At,  Bt]  Ap,  Bp\  and  Api,Bpi).  The  following  two  ionic  displacement 
mode  patterns  emerge  for  the  hybrid  in  Eq.(33)  taking  the  Bloch  superlattice  wavevector 
qSL  =  0  and  the  approximation  tanh.(q0dZns)  ~  1.  Both  ionic  displacement  patterns  are 
confined  within  the  ZnS  layer,  -dZns/ 2  <  2  <  dZns/ 2-  For  the  first  type, 


Zlx  — 


2iBetqxXqx{sin(qLz) 

-[1  -  Pi  tanh{qxdZnS/2)}  sm{qLdZnS/2) , 
sinh^z) 


sinh(g0dZn5/2) 


-Pism(qLdZnS/2)- 


r>> 


cosh(qxdZns/2)' 

2Bet9xXqL{cos(qLz) 

-^-[1  -  Pi  tanh(,,dz„s/2)]  Ma.dzns/2)^:i^Z), 

cosh  (qxz) 


(46) 


QlQo 

pi  sm{qLdZnS/2) 
Ql 


smh(q0dZnS/2) 


cosh(qxdZns/2) 
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and  for  the  second  type, 

ux  —  2BeiqxXqx{cos{qLz) 


where 


uz 


coshfg  z) 

-[1  -  p,  coth(gIrfz„s/2)]  cos(qLdzns /%)  5inh(9odz°  s/2) 

f  .  /oN  cosh(gIz)  , 
p2  cos(qLdZnS/2) sinh{qxdzns/2)  h 

2iBetqxXqiJ{sm(qiz) 

coa*hd**m  cos 

+  —  P2  cos(qLdzns/2)  .  ,  J?I  \o\}’ 


smh(qxdzns/ty' 


Pi  = 


Pi  = 


cosh(<fedsi/2)  cosh(qxdznS  fo) 
d 

smh(qxdSi/2)  smh(qxdZns/2) 


d  =r  sinh(qxdsi/2)  cosh(qxdZns/2)  +  sinh(&;dzns/2)  cosh(^xdsf /2) , 
€pi 


ep2 

and  epl  and  ep2  are  the  permittivities  in  Si  and  ZnS  layers,  respectively,  with 

_  k>2  ~  Ulo 

€p2~e0Ou2-LJ^o' 


(47) 


(48) 


(49) 


To  illustrate  the  patterns  of  ionic  displacements  in  the  ZnS  layers  given  in  Eqs.(46) 
and  (47),  we  need  to  first  determine  values  for  qx,  qL,  and  q0.  To  do  so,  we  will  follow 
the  numerical  procedure  described  in  Section  III  by  arbitrarily  fixing  a  value  for  the  in¬ 
plane  phonon  wavevector  qx  =  ir/(5azns),  where  aZns  is  the  lattice  constant  of  ZnS.  In 
calculating  the  carrier-optical  phonon  interaction,  the  value  of  qx  is  actually  determined 
by  the  conservation  of  in-plane  momentum  between  the  initial  and  final  states  of  the 
scattering  process.  For  a  given  value  of  qx,  typically,  a  set  of  hybridized  modes  can  be 
obtained.  Here,  we  show  only  the  mode  pattern  with  frequency  close  to  ojlo- 

We  obtained  fuo  —  41meV,  qi,  =  0.46  x  10 8/cm  and  q0  =  0.48  x  10s /cm.  Substi¬ 
tuting  these  values  into  Eqs.(46)  and  (47),  we  obtained  Figs.7(a)  and  7(b)  showing  the 
mode  patterns  of  ionic  displacement  of  both  the  first  and  second  types,  respectively  in 
a  ZnS  layer  of  dzns  =  20A.  It  can  be  seen  from  Figs. 7(a)  and  7(b)  that  the  mechanical 
boundary  condition,  vanishing  of  the  ionic  displacements  at  the  interfaces  of  Si  and  ZnS 
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layers,  is  satisfied. 


Figure  7.  Vibration  patterns  in  a  ZnS  layer  with  a  width  of  20 A  for  (a)  the  first  type 
and  (b)  the  second  type  solutions  of  the  hybridized  LO,  TO  and  IP  modes. 

III.  Dispersion  Relationship 


The  phonon  frequency  in  the  ZnS  layers  is  determined  by  the  following  set  of 
equations; 

f  u  =  uLO-  PM  +  ql), 

<  W  =  uto  -  Pr(ql  -  ql),  (50) 

„  h  +  h  cos (qLdzns)  +  h  sin (qL,dZns)  =  0 

where 

h  =  4psmh.(qxdSi)  +  4prsmh(qxdZnS), 

*2=  -4 pa, 

h  =  Sp2rsmh(qxdzns)sinh(qxdsi)  -  4p2a2  ^  ' 

+4 p2r2  sinh  2{qxdZns)  +  4 p2  sinh  2\qxdSi)  +  1, 


with 

and 


^  Aqirsd 

.  .2  ,  ,2 

UJ 

u LO  ~  uTO 


(52) 

(53) 


The  third  equation  in  (50)  is  obtained  from  the  requirement  of  a  nonzero  solution  for 
the  eight  simultaneous  equations  discussed  above,  and  Eq.(51)  is  arrived  at  under  the 
approximation,  tanh(g0cfens)  «  1. 


The  numerical  procedure  for  determining  a  phonon  frequency  is  the  following: 
given  a  value  of  qx,  we  can  determine  those  of  ty,  t2,  and  t3  from  Eq.(51).  Then  uj  is 
scanned  from  u>to  to  u>lo-  For  a  given  value  of  u>,  qL  and  qa  are  obtained  from  the  first 
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two  equations  in  (50).  Those  values  are  then  substituted  into  the  third  equation  in  (50) 
to  determine  if  the  particular  value  of  a;  is  a  solution. 

IV.  Scalar  and  Vector  Potentials 


Associated  with  the  two  types  of  ionic  displacement  in  Eqs.(46)  and  (47),  the 
scalar  potentials  in  ZnS  layers  are  given  as,  for  the  first  type, 


2p0BeiqxXsiu{qLZi)  \zx  \  <  ~~  ZnS  layer, 


I  *2 1  <  ^7jT  Si  layer, 


and  for  the  second  type, 


-2ip0BeiqxX cos(qLZi)  \zY\  <  ——  ZnS  layer, 

<j>  =  <  j  *  (55) 

0  1 22 1  <  -ir  Si  layer, 

k  z 

.Note  that  we  have  used  two  different  coordinates  zy  and  z2  for  layers  ZnS  and  Si,  re¬ 
spectively,  with  their  origins  placed  at  the  centers  of  the  respective  layers. 

The  vector  potentials  can  be  obtained,  for  the  first  type, 

mfeilM/i)  W<%£  ZnS  layer, 

Ax  =  ^  cosh{qxdzns/2)  2  (56) 

^MlBeiq*xVy  sinh($xz2)  \z»\  <  ^  Si  layer, 

^  cj  Z 

sinfefe^/2)  C°f (feZl)  N  <  %*  ZnS  layer, 

Az  =  <  .u  cosh{qxdzns/2)  2  (57) 

^^BeiqxXVy  cosh(qxz2)  |z2|  <  %  Si  layer, 

u  2 

and  for  the  second  type, 

'  cosfai  W2)  W<%*  ZnS  layer, 

^4  =  J  w  smh(9xdz„s/2J  2 

—  Be*gjXy2  cosh(gxz2)  |z2|  <  ^  Si  layer, 


4,= 


•  ^M^BeiqxXp2  cos{qLdZns/2)  .  ^ ^  <  ^  ZnS  layer’ 

a;  sinh(9xa^n5/2j  2 


4?xP0 


BeiQxXV2  sinh(7xz2) 


1*2  |  <  ~  Si  layer, 
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where 


Vi 


V2 


rsin(g£dZrt5/2)  cosh(gxdZnS/2) 

r  cosjgtdznsl^)  sinh(gadZng/2) 
2d 


and  pi,  P2,  and  d  are  given  in  Eq.(48). 


(60) 


The  scalar  potentials  associated  with  the  LO  modes  are  strictly  confined  within 
the  ZnS  layers.  Their  distributions  are  shown  in  Fig.8  for  the  first  and  second  types 
given  in  Eqs.(54)  and  (55)  with  qL  =  0.31  x  108/cm,  dzns  —  20 A,  respectively. 

The  vector  potential  associated  with  the  IP  modes  are  distributed  in  both  Si  and 
ZnS  layers,  even  though  Si  layers  are  treated  as  infinitely  rigid  and  do  not  contain  ZnS 
ionic  displacements.  The  profiles  for  the  two  components  of  the  vector  potentials  given 
in  Eqs.  (56-59)  for  the  first  and  second  types  with  dsi  =  40 A,  dzns  =  20A  are  shown  in 
Figs.9(a)  and  9(b),  respectively. 


Figure  8.  Scalar  potential  distribution  associated  with  the  LO  modes  in  a  period  of 
the  Si/ZnS  superlattice  with  dsi  =  40 A  and  dzns  =  20 A  for  both  the  first  and  second 
types  of  the  vibration  modes. 

It  can  be  seen  from  Figs. 8  and  9  that  both  scalar  and  vector  potentials  are 
not  continuous  across  the  interfaces.  However,  as  pointed  by  Ridley [9],  the  energy  of 
interaction  with  an  electron  traveling  coherently  with  the  optical  phonon  is  continuous. 
The  electric  field  can  be  obtained  as 

„  dA 

E  =  -V<f>- 

The  continuity  of  Ex  and  Dz  —  e(u>)Ez  implies  that 

wAx\z2=.±dsil2  ~  ~Qx<f>\zi=^dZns /2  d 
Az  |z2=±d5v/2  ^Az\Zl—~dz„s/^’ 


at  the  boundaries, 
-  wAx\Zl-zzdZnSl2) 


(61) 


(62) 
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where  A\x  and  A\z  are  x-  and  z-components  of  the  vector  potential  in  Si  layers.  The 
interaction  in  the  Si  layer  is  e(Aixvx  + Alzvz)  and  in  the  ZnS  layer  e(-<f>  + Axvx  +  Azvz), 
which  are  equal  when  the  electron  velocity  vx  =  w/qx  and  vz  =  0.  Thus,  the  coherent 
interaction  energy  is  continuous  across  the  interfaces. 


Figure  9.  Vector  potentials  associated  with  the  IP  modes  distributed  in  a  period  of  the 
Si/ZnS  superlattice  with  dsi  =  40A  and  dzns  =  20 A  for  (a)  the  first  type  and  (b)  the 
second  type  of  the  vibration  modes. 


The  electric  field  distributions  for  Ex  and  e(w)  in  Si  (d$i  =  40A)  and  ZnS  (dzns  = 
40A)  layers  are  shown  in  Figs. 10(a)  and  10(b)  for  the  first  and  second  types,  respectively. 
The  continuity  of  Ex  and  Dz  across  the  Si  and  ZnS  interface  according  to  Eq.(62)  is 
clearly  demonstrated. 


Figure  10.  The  field  distributions,  Ex  and  Dz,  derived  from  the  scalar  and  vector 
potentials,  in  a  period  of  the  Si/ZnS  superlattice  with  dSi  =  40A  and  dzns  =  20A  for 
(a)  the  first  type  and  (b)  the  second  type  of  the  vibration  modes. 
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V.  Summary  and  Discussion 


We  have  provided  an  analytical  model  of  optical  modes  in  Si/ZnS  superlattices 
consisting  of  polar  and  nonpolar  optical  phonons.  This  is  a  new  procedure  for  obtaining 
the  eigen  modes  of  a  mixed  polar-nonpolar  heterosystem.  In  the  Si  layers,  a  continuum 
model  with  double  hybridization  of  the  LO  and  TO  modes  is  used  to  describe  the 
vibration  patterns.  Since  there  is  no  electric  field  resulting  from  the  nonpolar  ionic 
displacements  in  Si  layers,  the  only  boundary  condition  that  needs  to  be  satisfied  in  the 
Si  layers  is  the  vanishing  of  the  displacements  at  the  Si-ZnS  interface,  as  the  ZnS  layers 
can  be  considered  as  infinitely  rigid  with  respect  to  the  vibrations  of  the  Si  layer.  Due  to 
this  strict  confinement,  only  guided  modes  emerge  in  the  Si  layers  which  consist  of  s-TO 
and  coupled  p-TO  and  LO  modes,  with  no  interface  modes.  These  guided  modes  have 
been  illustrated.  Their  interaction  with  carriers  in  the  superlattice  can  be  calculated 
through  the  optical  deformation  potential  for  Si.  The  interaction  Hamiltonian  can  be 
obtained  by  talcing  the  product  of  this  potential  with  the  normalized  ionic  displacement. 

However,  for  the  optical  phonons  in  ZnS  layers,  we  need  to  include  the  electrical 
interaction  in  calculating  the  carrier  scattering  by  optical  phonons,  since  there  are  elec¬ 
tric  fields  associated  with  the  polar  optical  vibrations.  As  a  result,  both  mechanical  and 
electrostatic  boundary  conditions  need  to  be  satisfied  in  the  interfaces.  A  continuum 
model  employing  a  linear  combination  of  LO,  TO  and  IP  (interface  polariton)  modes 
with  a  common  frequency  is  used  to  describe  the  ionic  displacements  in  ZnS  layers.  A 
numerical  procedure  for  determining  a  phonon  frequency  is  provided.  This  hybridized 
model  is  necessary  to  meet  the  simultaneous  requirement  on  the  mechanical  and  electro¬ 
static  boundary  conditions.  The  mechanical  boundary  condition  is  again  the  vanishing 
of  the  optical  displacements  since  Si  layers  can  be  considered  as  infinitely  rigid  with 
respect  to  the  vibrations  of  the  ZnS  layers.  The  electrostatic  boundary  conditions  are 
the  continuity  of  the  electric  field  parallel  to  the  interface,  and  the  continuity  of  the 
displacement  field  normal  to  the  interface.  Based  on  this  set  of  boundary  conditions, 
expressions  are  obtained  for  the  ionic  displacements  in  ZnS  layers  consisting  of  LO,  TO, 
and  IP  modes.  There  are  scalar  and  vector  potentials  associated  with  the  LO  and  IP 
modes,  respectively,  but  no  electric  field  associated  with  the  TO  mode.  The  scalar  po¬ 
tential  and  its  associated  electric  field  due  to  the  LO  mode  are  distributed  only  within 
the  ZnS  layers  and  are  zero  in  the  Si  layers.  But  the  vector  potential  and  its  associ¬ 
ated  electric  field  due  to  the  IP  mode  have  distributions  in  both  ZnS  and  Si  layers  even 
though  there  is  no  ZnS  ionic  displacement  mode  in  the  Si  layers.  Examples  of  these 
mode  characteristics  have  been  demonstrated.  Neither  the  scalar  nor  vector  potential 
is  continuous  across  the  Si-ZnS  interface,  but  the  energy  of  coherent  interaction  with 
carriers  is  continuous  due  to  the  continuity  of  the  electric  field  parallel  to  the  interface. 
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